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Bimetallic nanoalloys such as nanoparticles and nanowires are attracting significant attention 
due to their vast potential applications such as in catalysis and nanoelectronics. Notably, Pd-Pt 
nanoparticles/nanowires are being widely recognized as catalysts and hydrogen sensors. Compared 
to unary systems, alloys present more structural complexity with various compositional configura¬ 
tions. Therefore, it is important to understand energetically preferred atomic structures of bimetallic 
nanoalloys. In this study, we performed a series of simulated annealing Monte Carlo simulations 
to predict the energetically stable atomic arrangement of Pd-Pt nanoparticles and nanowires as a 
function of composition based on a set of carefully designed empirical potential models. Both the 
Pd-Pt nanoparticles and nanowires exhibit quasi-ordered configurations, quite similar to bulk alloy 
phases such as the Llo and the LI 2 structures with minor surface segregation effects. We believe 
that this study can provide a theoretical guide for the design of various bimetallic nanomaterials. 

PACS numbers: 61.46.-w, 64.70.Nd, 81.07.-b, 07.05.Tp 


I. INTRODUCTION 


II. COMPUTATIONAL DETAILS 


Recently, nanomaterials are drawing increasingly more 
attention due to their unique physical and chemical prop¬ 
erties that are distinct from bulk materials. Among var¬ 
ious nano-structured materials, bimetallic nanoparticles 
and nanowires are receiving the spotlight in many ar¬ 
eas such as catalysts, sensors, electronic devices and bio- 
medicals^ In particular, Pd-Pt bimetallic nanoalloys are 
well known as hydrogen sensors and catalysts—^ 

Numerous efforts have focused on the synthesis of 
bimetallic nanoparticles/nanowiresjA-. but their pre¬ 
ferred equilibrium structures are not yet well understood. 
Since the chemical and physical properties of bimetallic 
nanowires are mainly dependent on their size, composi¬ 
tional configuration and structural stability, it is impor¬ 
tant to identify their stable atomic structures. 

Compared to unary systems, nanoalloys present more 
structural complexity because the two components can 
exhibit various structural modifications. While bimetal¬ 
lic nanoalloys of various structures can be synthesized 
by a wide variety of techniques, the intrinsic equilibrium 
structure of bimetallic nanoparticles/nanowires depends 
on the alloy components. In spite of exhaustive works 
of both experiment and theory, it is not well established 
which alloy systems tend to exhibit mixed-alloy or or¬ 
dered intermetallic nanoalloys. While their structural 
and phenomenal complexity makes it difficult to experi¬ 
mentally predict the stable atomic structure of bimetal¬ 
lic nanoalloys ^ atomistic simulation could be used as an 
alternative to understand structural properties. In this 
paper, we performed a series of simulated annealing (SA) 
Monte Carlo (MC) simulations based on a set of embed¬ 
ded atom method (EAM) models to investigate energet¬ 
ically preferred atomic structures of Pd-Pt nanoparticles 
and nanowires. 


A. Interatomic potentials for Pd-Pt alloys 


One of the key components of atomic scale simulations 
is the method of evaluating the potential energy of sys¬ 
tems and interatomic forces as a function of the positions 
of atoms. Here, we adopted a set of original and updated 
versions of the embedded atom method (EAM) potential 
models for Pd-Pt developed in our previous work— 

The total energy of the system is given by the usual 
EAM form: 


*= E 
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where F Si (p) is the energy associated with embedding 
atom of type s, in a uniform electron gas of density p 
and <f> Si - Sj (r) is a pairwise interaction between atoms of 
type Si and s ? separated by a distance r. The electron 
density is given by 


Pi = X ( 2 ) 

j(Ai) 

where the atomic electron density function f(r) is taken 
as the density of a hydrogenic 4s orbital following the 
formalism by A. Voter—: 


f{r) = fo r 6 (e _/3r + 2 9 e" 2 ^) . (3) 

Here, (3 is an adjustable fitting parameter that quantifies 
the distance over which the electron density decays away 
from an atom position and fo is a prefactor. (We chose 
fo to be 1/ f{r eq Ni st ) for convenience.) 
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TABLE I: Parameters for the Pd-Pt EAM potentials. 



P (A- 1 ) 

Dm (eV) 

«M (A 4 ) 

VM (A) 

<5 (eV) 

Pd 

3.400 

1.682 

1.501 

2.343 

0.06 

Pt 

3.580 

0.783 

1.779 

2.577 

0.01 

Pd-Pt (I) 

- 

1.259 

1.701 

2.448 

0.00 

Pd-Pt (II) 

- 

1.254 

1.696 

2.443 

0.06 


The pair potential term <f>(r) is chosen to take a Morse 
potential form with a minor additional term: 


<t>{r) = 


—D m 
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+ (^? (r_rj) <r3 ‘ r) 
x0(r - r 2 )0(r 3 - r), (4) 


where Dm, Q-m, I'm, and S are adjustable fitting param¬ 
eters. Here, the second term in the pair potential func¬ 
tion was introduced in order to tune the melting point 
for unary metals and the heat of mixing for alloys: the 
parameter <5 represents the magnitude of pair interaction 
tuning between the second and third nearest neighbor 
positions, r 2 and r 3 , with the aid of the Heaviside step 
function 9(r). 

Finally, the embedding function was numerically deter¬ 
mined so that the total energy of reference unary system 
as a function of dilation satisfies the following universal 
binding energy relation^: 
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with a = \/9 BD/Eq, where a is the dilated lattice con¬ 
stant, a 0 is the equilibrium lattice constant and Eq, B , 
and H are the cohesive energy, bulk modulus, and equi¬ 
librium atomic volume of the reference lattice, respec¬ 
tively. The potential interactions were smoothly cut off 
at r = r cu t (~ 5.5 A) to ensure that the interatomic 
potential and its first derivatives are continuous. 

For a binary alloy potential, we need to fit the five ad¬ 
justable parameters (Dm, Rm, &m, P, and S) for each sin¬ 
gle component and the four parameters (Dm, Rm,&m, 
and 5) for a cross-species pair interaction. We adopted 
the parameter sets for Pd and Pt developed in our previ¬ 
ous work^ but re-optinrized the parameters for the Pd-Pt 
cross-species interaction by including the heat of mixing 
behavior of the solid and liquid phases and phase dia¬ 
gram information in the target properties. Table U shows 
the parameter set for the Pd-Pt system including two dif¬ 
ferent versions of Pd-Pt interactions, which were labelled 
I and II and hereafter referred as EAM-I and EAM-II, 
respectively. 



Pd 


Mole fraction of Pt 


Pt 


FIG. 1: Heat of formation for the Pd-Pt bulk system as 
a function of Pt concentration: The solid (black) and dashed 
(red) lines are the results of the present EAM models for a dis¬ 
ordered solid solution (random mixing) while the open circles 
are those of experimental measurement (Ref. EH). The for¬ 
mation energies of three possible intermetallics (LD-PdsPt, 
Llo-PdPt, and LH-PdPta) are also represented by solid sym¬ 
bols: The squares and diamonds with dotted lines correspond 
to the results of the present EAM models and the other solid 
symbols to first-principle calculation data from Refs. EH and 


B. Structural properties of Pd-Pt bulk alloys 

Compared to unary systems, alloy systems present 
more structural complexity due to various compositional 
configurations. Since atomic structures of nanomaterials 
are generally determined by competition between bulk 
mixing behavior and surface effects, the bulk alloy struc¬ 
ture and surface energy seem to be two key factors to 
be considered. However, in spite of being an apparently 
simple binary system, only a few old articles discuss its 
thermodynamic properties and bulk alloy phase behavior 
of the Pd-Pt system . 11 : 12 

Figure Q] shows the calculated heats of formation for a 
disordered solid solution (random mixing) and some pos¬ 
sible intermetallics in comparison with experimental and 
first-principle calculation data on the Pd-Pt system in 
the literature ; 11 : 13 : 14 At higher temperatures, the Pd-Pt 
system forms a completely miscible solid solution based 
on the face-centered cubic (fee) lattice^ and both the 
EAM-I and EAM-II potentials agree well with the exper¬ 
imental data. However, the Pd-Pt system has not been 
studied in detail at low temperatures, and there has been 
controversy over whether or not any stable intermetallic 
compounds exist. For example, most experiments and 
some first-principle calculations based on local density 
approximation (LDA) propose that Llo is the most sta¬ 
ble compound for Pd-Pt of 1:1 stoichiometry , 13 i 14 while 
some other LDA calculations predict that Lli structure 
has the lowest energy ; 15 : 16 

Table Ell provides the calculated formation energies of 
some fee based alloy crystal structures in comparison 
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TABLE II: Calculated enthalpies of formation for the pos¬ 
sible Pdi-zPtz intermetallics: the formation energies of 
some face-centered cubic (fee) based alloy crystal structures 
(Llo, Lli, I/I 2 , DO 22 ) are given in comparison with those by 
first-principle calculations. The numbers with * indicate the 
formation energies of the lowest energy structure. 


Crystal Structure 

EAM-I 

EAM-II 

lda±£ 

LDAf^ 

Llo {PdPt) 

-68.7* 

-64.2* 

-39.5* 

-77.7* 

Lli {PdPt) 

-24.9 

-26.3 

-29.1 

-68.4 

LI 2 {Pd 3 Pt) 

-50.0 

-46.8* 

-30.4* 

-58.8* 

V0 22 {Pd 3 Pt) 

-50.6* 

-44.0 

-22.1 

-54.4 

LI 2 {PdPt 3 ) 

-48.2 

-45.1* 

-35.1* 

-57.8* 

m 22 (PdPt 3 ) 

-48.7* 

-41.9 

-29.1 

-51.3 
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FIG. 2: Calculated phase diagram of the Pd-Pt binary al¬ 
loy (solid lines) in comparison with experimental data (solid 
symbols, Ref. fl2l) . 


with those by first-principle calculations. 13,14 ' As shown 
in Fig. [I] both the EAM potentials give reasonable val¬ 
ues for the formation energy of Llo and LI 2 intermetallic 
compounds falling in the middle of the range of Refs. [H 
and [lj. As indicated in Table El the formation energies 
of the DO 22 (PdsPt and PdPta) structures are lower than 
that of LI 2 structure (by less than 1 meV) in the case of 
the EAM-I potential, while the EAM-II model predicts 
the (100) family of ordered states that consists of LI 2 for 
PdaPt and PdPt 3 and Llo for PdPt in accordance with 
Refs. [3 andflil 

We also calculated the solid-liquid phase diagrams 
for the Pd-Pt bulk alloys using Gibbs-Duhem integra¬ 
tion technique combined with semigrand canonical Monte 
Carlo simulations! 17 ’ 18 Figure [2] shows calculated phase 
diagrams of the Pd-Pt binary alloy in comparison with 
experimental data.— Although the parameter set of the 
EAM-II potential was marginally modified from that of 
the EAM-I potential, it produced results almost identical 
to those produced by the EAM-I potential in predicting 
the enthalpy of formation for a disordered solid solution 
(random mixing). However, their solid-liquid phase di¬ 
agrams are quite different: Since the potential parame¬ 



FIG. 3: Initial atomic configurations: bimetallic nanoparti¬ 
cles of random mixing with (a) truncated octahedron (TOh) 
and (b) icosahedron (Ih) shapes and bimetallic nanowires of 
random mixing with (c) [100], (d) [110], and (e) [111] crystal¬ 
lographic orientations along the nanowire axis. 


ters are adjusted to reproduce the heats of formation for 
both the solid and liquid phases, the solid-liquid phase 
diagram predicted by the EAM-II model agrees much 
better with the experiments. 

In nearly all of the previous simulation studies on Pd- 
Pt nanoalloysa simplified combination of the Pt-Pt 
and Pd-Pd parameters was set as potential parameters 
for the cross-species interaction instead of fitting to Pd- 
Pt alloy properties. Consequently, bulk Pt-Pd alloys were 
assumed to be solid solutions, rather than ordered inter- 
metallics, for all compositions.—^ 6 However, our EAM 
models for Pd-Pt alloys, which were fitted to alloy prop¬ 
erties such as the formation enthalpy and the lattice con¬ 
stant, rather predict stable intermetallic compounds in 
accordance with the first-principle calculations. 

Although the bonding characteristics of materials orig¬ 
inate from their electronic structures, in classical atom¬ 
istic simulations, such as molecular dynamics (MD) and 
Monte Carlo methods, the interactions between atoms 
are usually described by empirical potentials, and there¬ 
fore the reliability of the simulation results is entirely 
dependent on the reality and accuracy of the interatomic 
potentials. So far, lack of appropriate potential models 
for bimetallic systems have attributed to insufficiencies 
of atomic scale simulations on the structural properties 
of alloy nanoalloys. Considering that the present EAM 
potentials for the Pd-Pt alloy system well reproduce the 
formation energies of the alloy phases as well as the ex¬ 
perimental phase diagram, the present empirical model 
approach may provide reliable simulations of alloy nanos¬ 
tructures such as nanoparticles and nanowires. 


C. Monte Carlo simulations of nanoalloys 

It is well known that metallic nanoparticles present 
competitive structural motifs, such as the icosahe- 
dron(Ih), decahedron(Dh), cuboctahedron, and trun¬ 
cated octahedron (TOh). 27 ' 28 In our atomistic simula¬ 
tions of alloy nanoparticles, truncated octahedral (TOh) 
and icosaliedral (Ih) nanoparticle of ~ 4 nm in diam¬ 
eter (1654 and 1415 atoms, respectively) were used as 
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initial atomic configurations. Initial compositional con¬ 
figurations were typically generated by random mixing of 
atoms as shown in Figs. HHa) and (b). 

For the simulations of nanowires, cylindrical shape 
nanowires of ~ 4 nm in diameter (also with random 
mixing) were used as initial structures. The bimetallic 
nanowires were chosen to have single-crystalline fee struc¬ 
tures with three different crystallographic orientations, 
i.e., [100], [110], and [111] along the nanowire axis, as 
shown in Figs. Hc)-(e). We also investigated the bimetal¬ 
lic nanowire structure of various mole fractions: the com¬ 
position of 25, 50 and 75 at% Pt. The nanowires were 
modelled as infinitely long wires by applying the periodic 
boundary condition along the wire axis. 

In order to investigate energetically preferred atomic 
structure of Pd-Pt nanoalloys, we carried out a series of 
simulated annealing Monte-Carlo optimizations based on 
Metropolis algorithm^ The simulated annealing proce¬ 
dure was performed from T = 1000 K to 0 K at a cooling 
rate of ~ 1 K/MCS (where MCS stands for Monte Carlo 
step) in order to avoid being isolated in a local minimum. 
In all of the MC simulations, lattice relaxation using the 
conjugate-gradient method was implemented after con¬ 
figurational exchange of two select atoms. 


III. RESULTS AND DISCUSSION 

A. Structure of Pd-Pt bimetallic nanoparticles 

We found that the optimized structures were relatively 
insensitive to the initial structure or cooling rate of the 
simulated annealing procedure, although there were some 
minor changes in surface atomic structure due to the ran¬ 
dom nature of Monte Carlo simulations. The consistent 
results may be attributed to the low transition barrier 
for compositional configuration change from a random 
structure to the optimized structure in the simulations. 

Figure [4] shows the typical optimized atomic struc¬ 
tures of Pd-Pt nanoparticles of 1:1 stoichiometry that 
were simulated with the EAM-I and the EAM-II poten¬ 
tials. Both the EAM potential models produced similar 
inner structures with some minor difference in surface 
atomic structure: The cross-sectional views of truncated 
octahedron Pd-Pt nanoparticle depict a typical Llo inter- 
metallic compound structure with surface segregation of 
Pd, while alloy nanoparticles of icosahedral shape exhibit 
somewhat complicated patterns in their cross-sectional 
views, probably due to the intrinsic twin boundaries of 
the icosahedral nanoparticles. Note that all of the {100} 
type facets of TOh nanoparticles are occupied only by 
Pd atoms, while {111} type facets are composed of both 
Pd and Pt atoms with linearly aligned patterns. In both 
the TOh and Ih nanoparticles, edges are mostly occupied 
by Pd atoms. 

In order to clearly understand the inner structures, the 
bond characteristics of the alloy nanoparticles were ana¬ 
lyzed by using the short-range-order parameter proposed 


(a) 






FIG. 4: Predicted atomic structures of Pd-Pt nanoparticles 
simulated with (a) the EAM-I and (b) the EAM-II potentials: 
The first and second columns of the figures show each of their 
surface and cross-sectional views of the truncated octahedron 
(TOh) nanoparticles while the third and fourth columns of the 
figures show each of their surface and cross-sectional views of 
the icosahedron (Ih) nanoparticles, respectively. Pd atoms 
are in blue (dark) and Pt in grey (light). 
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FIG. 5: Short-range order parameter ai for the TOh and Ih 
nanoparticles in comparison with those of the ideal ordered 
alloy structures. The order parameter at was calculated for 
the first three shells of nearest-neighbors(NN). 


by Cowley^H and Warren£L The Warren-Cowley short- 
range order parameter on is defined as 

on = 1 ~Pi,AB/C B = 1 - (6) 

where Pi,AB{=Ni : AB/Ni) is the probability that the 
atomic site is occupied by a B atom in the ith-neighbor 
shell of an A atom, Ni^AB is the expected number of 
B atoms around an A atom in the ith shell of neighbors, 
Ni is the total neighbor number in the ith-neighbor shell, 
and C’B is the atomic concentration of B. All of the ct; 
values range between -1 to 1 and positive values denote 
more A-A and/or B-B pairs and negative values denote 
more A-B pairs between the central and the neighbors in 
the ith shell, compared to those of random mixing. The 
parameter for the first-nearest neighbor shell, on, is of¬ 
ten referred as chemical short-range order parameter and 
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is a useful measure of the chemical affinity representing 
the degree of tendency for ordering: positive aq indicates 
clustering or phase separation and negative oq indicates 
strong A-B bond and chemical ordering. 

Figure [5] shows the first three short-range order param¬ 
eters on for the TOh and Ih nanoparticles. (Here, cq were 
averaged over several configurations obtained in simula¬ 
tions performed with different random number seeds.) 
The oti for the ideal Llo and Lli ordered alloy structures 
were also given for comparison: For example, for the ideal 
Ll! lattice, a 2 = —1, while ctq and 0:3 are zero. For the 
Llo alloy structure, a-i = 1, while aq and 03 are 1/3. 
On the other hand, in the case of randomly mixed alloys 
(not shown here), oq tend to zero for any *-th neighbor 
shells. 

As indicated by the distributions of aq shown in Fig. [5] 
atomic structures of the TOh Pd-Pt nanoparticles simu¬ 
lated with both the EAM potentials exhibit very similar 
bond characteristics to that of the Llo alloys: both aq 
and 0:3 are negative toward -1/3 and 0:2 is strongly posi¬ 
tive. The distributions of at for the Ih nanoparticles also 
show similar patterns as that of the TOh nanoparticles 
(indicating multiply-twinned Llo structures), although 
the absolute value of ai is relatively small. Interestingly, 
while the values of a\ are almost identical in both the 
EAM models, 012 and 0:3 by the EAM-II model are closer 
to 1/3, the ideal value for the Llo structure, which in¬ 
dicates that the nanoparticle structures predicted by the 
EAM-II model are more ordered with longer range-order. 


B. Structure of Pd-Pt bimetallic nanowires 

In the case of the Pd-Pt nanowires, both the EAM-I 
and EAM-II potentials predicted nearly identical views 
of the structures. Therefore, we limit our focus to the 
simulations based on the EAM-II potentials (of more ac¬ 
curate phase diagram behavior) in the remainder of this 
paper. 

Figure [ 6 ] shows the typical optimized atomic structures 
of Pd-Pt nanowires with three different crystallographic 
orientations, i.e., [ 100 ], [ 110 ] and [ 111 ] along the nanowire 
axis. In all three cases, the inner patterns of alternating 
distribution of Pd and Pt atomic planes appear to be a 
minor modification of the Ll 0 alloy structure with Pd 
atoms segregated on the surface. Interestingly, the alter¬ 
nating {100} type planes of Pd and Pt atoms are mostly 
aligned parallel to the nanowire axis in the case of [ 100 ] 
oriented nanowires probably because it reduces the sur¬ 
face energy and internal stresses. In the [110] and [111] 
oriented nanowires, alternating (100) planes of Pd and Pt 
are inclined according to their crystal orientation angles. 

FiguresjTJa) and (b) show the typical optimized atomic 
structures of Pd-Pt nanowires with compositions of 25 
and 75% of Pt, respectively. In the case of 3:1 or 1:3 
compositions, the internal structures of nanowires show 
that minor element atoms occupied the lattice sites with 



FIG. 6 : Predicted atomic structures of Pd-Pt nanowires with 
(a) [ 100 ], (b) [ 110 ], and (c) [ 111 ] crystallographic orientations 
along the nanowire axis. The first column of the figures shows 
the tilt view of the nanowires, while the second, third, and 
last columns of the figures show cross-sectional views of the 
xz-, yz-, and xy- planes, respectively. Pd atoms are in blue 
(dark) and Pt in grey (light). 



FIG. 7: Predicted atomic structures of Pd-Pt nanoparticles 
with compositions of (a) 25%, (b) 75% Pt. The first col¬ 
umn of the figures shows the tilt view of the nanowires, while 
the second, third, and last columns of the figures show cross- 
sectional views of the xz-, yz-, and xy- planes, respectively. 
Pd atoms are in blue (dark) and Pt in grey (light). 


regular patterns. Although the minor element sites do 
not show perfect long-range order, quasi-ordered atomic 
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FIG. 8: Short-range order parameter on for the nanowires in 
comparison with those of the ideal ordered alloy structures. 
The order parameter on was calculated for the first three shells 
of nearest-neighbors(NN). 


arrangement appears to be that of the LI 2 bulk alloy 
phases, as shown in the cross-sectional views of Fig. Q] 
Compared to the distribution of Pt atoms on the surface 
of 25% Pt nanowires, Pd atoms are strongly segregated 
to the surface for the nanowires of 75% Pt presumably 
due to the lower surface energy of Pd. On the other 
hand, Pt atoms tend to be segregated in the layers be¬ 
neath the surface (instead of the surface layers) at its low 
concentration, as shown in Fig. [7[a) 

Figure [ 8 ] shows the short-range order parameters on 
calculated for the optimized structures of the nanowires. 
The on for the ideal Ll 2 (same as Ll 0 ) and DO 22 ordered 
alloy structures were also given for comparison. As in¬ 
dicated by the distributions of oti shown in Fig. [ 8 ] all of 
the Pd-Pt nanowires of 50% Pt exhibit very similar bond 
characteristics to that of the Llo alloys with ai <-0.25, 
«3 <-0.2, and a 2 >0.2. The distributions of at for the 
25% and 75% Pt nanowires also show similar patterns as 
those of the ordered structures, although it cannot be ab¬ 
solutely identified as either the Ll 2 or the DO 22 structure. 
Considering that their inner lattice does not satisfy 3:1 or 
1:3 compositions due to surface segregation effects, their 
internal structures are more like those of ordered struc¬ 
tures at various diameters up to several nanometers. 

Generally, equilibrium compositional structure of 
bimetallic nanoalloys may be either similar to the bulk 
alloys or different from the bulk by the competition be¬ 
tween atomic mixing behavior and surface effects. Ta¬ 
ble mu shows the calculated surface energies for various 
surface types in the Pd-Pt system. Although the (100) 
surface energy of Pd is substantially lower than that of 
Pt, the (111) surface energy of Pd is comparable to that 
of Pt. The surface energy values of randomly-mixed al¬ 
loys are between those of pure Pd and Pt and the surface 
energies of the Llo alloy are even lower than those of dis¬ 
ordered alloys. Because of the comparable surface energy 


TABLE III: Calculated surface energies for various surface 
types in the Pd-Pt system. 


Bulk structure 

Surface 

(100) 

energy [mj/m 2 ] 

(111) 

Pd (fee) 

1815 

1700 

Pt (fee) 

1891 

1676 

PdPt (disordered) 

1847 

1678 

PdPt (Llo) 

1755 

1590 


properties between those surfaces, the Pd-Pt nanoalloys 
may maintain a similar compositional structure to the 
bulk ordered phases without dramatic surface modifica¬ 
tion, even at this nanoscale regime. 


IV. SUMMARY 

In experiments, energetics is not the only factor that 
affects the morphology of nanoalloys such as nanoparti¬ 
cles and nanowires. In cases where the formation is de¬ 
termined by kinetics rather than thermodynamic factors, 
the details of individual synthesis processes can govern 
the structures of the synthesized nanoalloys. Neverthe¬ 
less, identifying energetically stable compositional con¬ 
figurations may be a first step towards understanding 
structural stability. 

For detailed structural information of the energetically 
stable bimetallic nanoparticles and nanowires, we per¬ 
formed a series of Monte Carlo optimizations with sim¬ 
ulated annealing procedure based on a set of carefully 
designed EAM potential models. The EAM-I and EAM- 
II models predict fairly similar structures for nanoalloys, 
although they are distinguishable in their solid-liquid 
phase diagram behavior. Our simulations also investi¬ 
gated the energetically favored atomic structures of Pd- 
Pt nanoparticles and nanowires with various shapes and 
compositions. Both Pd-Pt nanoparticles and nanowires 
exhibit quasi-ordered configurations, quite similar to 
bulk alloy phases such as the Llo and the Ll 2 structures 
with minor surface segregation effects. 

Our simulations quantify the energetically favorable 
atomic arrangements of the Pd-Pt nanoalloys and explain 
why they exhibit such equilibrium structures in terms of 
their bulk alloy structures and surface energies. We be¬ 
lieve that our simulation results may be useful reference 
for further investigation of Pd-Pt nanoalloys. Above all, 
we expect that this study will provide a theoretical stan¬ 
dard for design of bimetallic nanoparticles and nanowires 
with various possible structures. 
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